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Recent progress in optical lattice clocks requires unprecedented precision in controlling systematic 
uncertainties at 10“^® level. Tuning of nonlinear light shifts is shown to reduce lattice-induced clock 
shift for wide range of lattice intensity. Based on theoretical multipolar, nonlinear, anharmonic and 
higher-order light shifts, we numerically demonstrate possible strategies for Sr, Yb, and Hg clocks 
to achieve lattice-induced systematic uncertainty below 1 x 10“^®. 

PACS numbers: 06.30.Ft, 32.60.-hi, 37.10.Jk, 42.62.Eh, 42.62.Fi 


INTRODUCTION 


Last few years have witnessed significant advances in 
optical clocks to reach uncertainties of 10“^® level in ion- 
based clocks [H and optical lattice clocks ii- Hitherto 
unexplored accuracy of optical clocks opens up new possi¬ 
bilities in science and technologies, such as probing new 
phwics via possible variation of fundamental constants 
[J-y, and relativistic geodesy to measure gravitational 
potential differences [ij . Evaluations of perturbations on 
the clock transitions are indeed at the heart of these en¬ 
deavors. 

Unperturbed transition frequencies have been accessed 
by extrapolating perturbations to zero, which is straight¬ 
forward if the correction is proportional to the perturber. 
Once the dependence becomes nonlinear, such as the 
blackbody radiation shift that changes as T'^ with tem¬ 
perature [^, dedicated experimental iii and theo¬ 
retical investigations 0 are crucial. Contrarily, non¬ 
linear response is leveraged to make the clock transition 
frequency insensitive to perturbations in hyper-Ramsey 
spectroscopy [Hi- 

Optical lattice clocks aimed at eliminating light-shift 
perturbations on the clock transition by operating an op¬ 
tical lattice at the magic frequency 12|, which at first 
glance should exempt them from evaluating lattice-laser 
intensity. However, residual light shifts aris ing from hy¬ 
perpolarizability [13| and multipolar effects IJ, Il5j man¬ 
ifest as leading systematic uncertainties at low 10“^^. 
Coupled with atomic motion in the optical lattice l4ll6| 
whose intensity varies in space to confine atoms, these 


light shifts show intricate nonlinear response to the lat¬ 
tice intensity 17|, which makes the corrections highly 
delicate issue. 

In this Letter, we propose strategies to eliminate light 
shift perturbations without relying on a conventional 
zero-extrapolation approach. We show that lattice laser 
frequency and feht-polarization-dependent hyperpolar¬ 
izability effect [1^ can be used to tailor intensity depen¬ 
dence of light shift. We define an “operational magic 
frequency” to reduce light shift to less than 1 x 10“^® for 
a sufficiently larger intensity variation than that is nec¬ 
essary for conhning atoms. Numerical calculations for 
electric-dipole (El), magnetic-dipole (Ml), and electric- 
quadrupole (E2) polarizabilities and hyperpolarizabilities 
are presented for the clock transitions in Sr, 

Yb, and Hg atoms, which are used to demonstrate the 
feasibility of the proposed strategies. 


LATTICE INDUCED LIGHT SHIFTS 

The lattice-induced light shift is given by the difference 
between the Stark energies in the ground (g) and excited 
(e) clock states. We assume a one-dimensional optical 
lattice as depicted in Fig. [T] and consider the interaction 
between a trapped atom and a standing-wave field, 

E(a;,t) = 2EoCOsfca;coswt, (1) 

which consists of counter-propagating laser waves with 
electric field vector Eq (corresponding intensity /), fre¬ 
quency w, and wavevector k = with k = uj/c and 
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FIG. 1. Clock transitions in an optical lattice potential cre¬ 
ated by counter-propagating laser beams. An atom trapped 
in Ug{x) with \x\ Am/4 is excited on the ^S'o(g) ^^’o(e) 

clock transition with unperturbed frequency vq and lattice- 
induced clock shift Ai/c. Actual laser beams show trans¬ 
verse intensity distribution with beam radius 

ro S> A, which weakly confines atoms transversely. 


c the speed of light. The atom-lattice interaction is de¬ 
termined by the operator V{x,t) = Re[y(x) exp(—iwt)] 
with the spatial factor 

V(x) = Vei cos kx + (Ve 2 + Vmi ) sin kx, (2) 

where Vei, Ve2, and 14ii correspond to operators of El, 
E2, and Ml interactions. 

The second- and fourth-order terms of atom-lattice in¬ 
teraction energy V(x) correspond to linear and quadratic 
terms in lattice-laser intensity I. The optical lattice po¬ 
tential for an atom at |a;| <g; A = 27r/k (see Fig. [T]) is 

given by [13 

with potential depth given by 

i?g(e)(/) = -C/g(e)(0,/) = + (4) 

which is determined by the electric-dipole polarizabil¬ 
ity and hyperpolarizability /3g(e){W)0 with ^ 

the degree of circular polarization of light as discussed 
later. The coefficient for the harmonic term in Eq. ©, 

determines the vibrational fre¬ 
quency f2g(e)(-^) of atoms in the lattice. In terms of the 
photon recoil energy = {hk)'^/2Ai with Ai the atomic 
mass and H = h/2TT the Planck constant, the vibrational 
frequency is given by [13 


where a combined E1-E2-M1 polarizability 




e) 


^g(e 


(5) 


( 6 ) 


is the difference between El polarizability and the sum of 
E2 and Ml polarizabilities ag™^(a;) = j(w)-|-a^gj(w). 


TABLE I. Numerical values of the electric-dipole (q;bi), dif¬ 
ference of multipolar (Aa'*™) and hyper (A/I) polarizabili¬ 
ties for Sr, Yb and Hg clock transitions at the magic wave¬ 
lengths Am. Frequency shifts and vibrational frequencies 
12(7) « f 7 are calculated for an atom trapped at 

the anti-node of the standing wave, where 7(kW/cm^) de¬ 
notes the intensity of each traveling wave laser. Non-zero 
imaginary part of the differential hyperpolarizability A/I*^'^^ 
for Hg accounts for two photon ionization rate in the optical 
lattice, where I and c correspond to linear and circular lattice- 
light polarization. Tr is the lattice photon recoil energy. A 
merit factor k indicates the insensitivity to atomic-motion- 
dependent lattice light shift. Operational lattice intensity 7op 
assumes to produce 5 times deeper lattice potential depth 
than the energy of laser-cooled atoms. 


Atom 

Sr 

Yb 

Hg 

Am (nm) 

813.4 

759.4 

362.6 

uo (THz) 

429 

518 

1129 


45.2 

40.5 

5.70 


1.38 

-1.71 

8.25 


-200 

-309 

-2.20 -1- 0.82i 


-311 

238 

4.40 -f 1.21i 

n f kHz j 

2n,/l y.^kW/cm2 J 

25.05 

18.03 

13.1 

1- 10^ ■ A 

hou V^kW/cm"^ / 

0.254 

0.720 

0.134 

Sn/h (kHz) 

3.47 

2.00 

7.57 

K = a®V|Aa'i”| 

3.3 X 10’’ 

2.4 X 10'^ 

6.9 X 10® 

7op = 5kBT/a^^ (kW/cm^) 

2.3 

10 

550 


The lowest-order anharmonic correction in the lattice is 
given by 


n%il)= <;‘,7(a;)/ + 5/3g(3)(.;,0/^ (7) 


dqm 


1 k^ 


J 3 


The energy of an atom in the n-th vibrational state |n) 
is calculated to be 


+ fi-flg(e) (I) (^ + l) ~^g{e) +’^+5)1 (8) 

where the second term corresponds to the harmonic- 
oscillator energy and the last term the anharmonic cor¬ 
rection given by 


p(4) f T\ _ ^R 


1 + 


3/3g(e)(^,C)7 

<reT(-) . 


(9) 


The lattice-induced clock light shift is given by the differ¬ 
ence of energies [Eq. (jS])] for the atom in its ground and 
excited states with vibrational state jn) being unchanged 
(the Lamb-Dicke regime). 


hAv,{I,n) = E:^'°{I,n) - E^^'°{I,n) = -AD{1) 

A hAn{I) (n + i) - AE;(4)( 7 ) {n^ + n+^), (10) 
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FIG. 2. The wavelength-dependent hyperpolarizability for Yb 
for linear (Ap'', dashed) and circular solid) polarized 

lattice near the magic wavelength Am. The vertical dashed- 
lines indicate two-photon resonances. 


where we define AD{I) = De{I) — Dg{I), An(/) = 
Qeil) - and - E^'^\l). 

The hyperpolarizabilities depend on the lattice-laser 
ellipticity in addition to its frequency, and are given by 




/3g(e)(w) -/3g(e)(w) 


( 11 ) 


where /3g(e) (w) is the hyperpolarizability for linear (cir¬ 
cular) polarized light. The degree of circular polariza¬ 
tion is defined by ^ = sin 2 %, where tan x determines 
the ratio of the minor to major axis of the polarization 
ellipse with ellipticity angle defined in 0 < y < 7 r/ 4 . 
When A/3* (= /3g — /3g) and A/3'^(= /3g — /3|) have opposite 
signs, there exists a “inagic ellipticity” determined by 
Cm = 1 / V^l — A/3'^/A/3* [I^ , which eliminates the differ¬ 
ential hyperpolarizability A/3(C) = A/3* — A/3*) 

[Eq. (fTTI) ]. However, a more important consequence for 
the following discussion is the tunability of A/3(C) be¬ 
tween A/3* and A/3°. 

Figure [2] shows the hyperpolarizability A/3*^'^^ for Yb 
atoms near the magic wavelength. The two-photon res¬ 
onance on the 6s6p^Po —>■ 6 s 8 p^Po transition at A 2 = 
759.71 nm appears only for linear polarized lattice, which 
locates between Ai = 754.23 nm and A 3 = 764.95 nm that 
correspond to two-photon resonances on the 6s6p^Po —>■ 
6 s 8 p^P 2 , 6 s 5 /^F 2 transitions [l^. These resonances al¬ 
low A/3* and A/3° to have opposite signs. Similar discus¬ 
sion holds for Hg, however, not for Sr. Table |T] summa¬ 
rizes the susceptibilities for Sr, Yb, and Hg clock tran¬ 
sitions calculated in the model-potential approximation 
[ 2 ^. While some of these values require experimental in¬ 
vestigations, in the following, we apply them to demon¬ 
strate the concept of the operational magic wavelength. 


OPERATIONAL MAGIC FREQUENCY 


In order to clarify the light-shift dependence on inten¬ 
sity I, we approximate the light shift [Eq. (fT0)l ] assum¬ 
ing experimentally feasible lattice laser intensities to trap 
laser-cooled atoms as listed in Table HI Following quanti¬ 
ties, (i) the E2-M1 polarizabilities -b 

(ii) the hyperpolarizability effect Pg(e)I, and (hi) the 
differential dipole polarizability Aa^^ = = 

— Suj, are about 10® times smaller than the electric- 
dipole polarizability <ag(e) ~ The light shift is then 
expanded in Taylor series in the vicinity of w ~ and 
neglecting higher order terms. 


hAvc{I, n, 5v, = Ci/ 2 . 1 ^^^ + ci/ -b -b 


dv 
SAa^i 


-5y - Aaa“ {2n -b 1) 




dv 


Sv + A/3(^) (2n^ -b 2n -b 1 ) 


4a®^ 
3fR 


jl/2 


4a^^ 


+ Ami‘2n+l)^I^I^/^-Ami\ ( 12 ) 


where 6v(= 5ui/2 t:) is detuning from the “El-magic fre¬ 
quency” defined by Aa^^(a;)(;^) = 0, Aa^™ = Q-q™ _ Q-q™ 
is the differential multipolar polarizability. 

The magic frequency Wni/27r(= c/Am) so far aimed 
at minimizing ci, which dominates Eq. (HU), by tuning 
Wm [l2|. However, this protocol is no longer valid 

for fractional uncertainty AvJvq ^ 10~^^ as the other 
CjP terms equally contribute [1^. We may define a merit 
factor K = a^^/|Aa‘i“| in Table U to indicate applica¬ 
bility of the “El-magic frequency” by neglecting multi¬ 
polar effect. K~^ indicates the fractional contribution 
of the E2-M1 effect that introduces atomic-motion in¬ 
duced nonlinearity. Since typical clock experiments 
ii are performed in a relatively low intensity regime 
(A/3J < Aa'i™), a large k for Sr and Yb may validate 
a linear extrapolation of the clock shifts down to 10“*^^ 
level. However, as the merit factor decreases for Hg, this 
simple approach breaks down and new strategies are re¬ 
quired. 

Figure [3] illustrates intensity dependence of the light 
shift Av^^{I) for Hg. For low laser intensity, behav¬ 
ior dominates, whose coefficient Ci /2 is determined by the 
electric-dipole and multipolar (Aa'^™) polar¬ 

izabilities. For intermediate intensity where clocks oper¬ 
ate, the leading term oc / is determined by the electric- 
dipole polarizability and slightly by the hyperpolarizabil¬ 
ity A/3(^) effect via the anharmonic correction [Eq. (|9])], 
in addition, the and P terms depend on A/3(^). 

It is apparent that there is no “magic condition” that 
sets four coefficients Cj 0, as there are only three 
free parameters n, Sv and available. Consider that 
an operational lattice intensity Pp is naturally deter- 
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FIG. 3. Light shift for a Hg clock as a function of laser inten¬ 
sity I. Black-solid and blue-dashed lines correspond to linear- 
polarized and “magic-elliptical” light with detuning 5v = 0. 
With 5v = —4.66 MHz and = 0.75, the light shift (red- 
solid) becomes less than ±1 mHz (green) for gray shaded re¬ 
gion. Red-dashed lines indicate tolerance (0.5%) for 5**®. For 
linear polarized light (^ = 0) with 5v = —2 MHz, the light 
shift becomes insensitive to A7 around I ~ 36 kW /err? (dot¬ 
ted line). 
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Lattice laser intensity / (kW/cm^) 

FIG. 4. Contour plots of light shifts for a) Yb and b) Sr 
clock transitions as functions of lattice laser intensity I and 
detuning 5v, for = 0.75 and = 0. The red-dotted 
lines show zero light shift. The region bound by red lines 
corresponds to light shift | Ar'c|/r'o < 1 x 10“^®. Tuning to the 
operational magic frequency 5v as indicated by white dotted 
lines, wide range of operational intensities are allowed. 


mined (Table H]) so as to trap Doppler-cooled atoms 
Dg(e) « a^^Iop ~ 5 with temperatures of T ~ 1 uK , 
4/iK, and 30 for Sr, Yb, and Hg, respectively (2ll - (^ . 
We may then tune Si' and ^ to minimize the light shift 
Ai/c for /op with its allowance A/ as wide as possible, 
as the laser intensity is spatially non-uniform. Actual 
one-dimensional optical lattices employ Gaussian beams 
with an intensity profileand a beam radius 
tq (3> A) to confine thermal motion of atoms in the trans¬ 
verse direction (see Fig. [I}, which inevitably introduces 
intensity inhomogeneity of A/// « 

In the following, to simplify the discussion, we assume 
atoms in the vibrational ground state (n = 0) along x- 
axis by applying a sideband cooling. We then look for 


conditions that satisfy both 


= 0 and 


I = Io 


Ar'c(/op, 5/^, C) = 0- The red solid line in Fig. |3] shows 
the light shift Ar'^s(/) with 5v = —4.66 MHz and = 
0.75, which demonstrates that the light shift becomes less 
than 1 mHz (corresponding to fractional uncertainty of 
1 X 10“^®) for 115kW/cm^ < / < 177kW/cm^. Here, 
the two photon ionization rate lva[?/h]I^p ~ 0.02 Hz for 
/op 150kW/cm^ is negligible for obtaining a few Hz 
linewidth. 

Figure!!^ calculates the light shift Ar'J’^(/, 5v) of Yb as 
functions of lattice laser intensity / and detuning Sv with 
„ 0.75, where a slight difference from the magic el- 
lipticity = 0.7516 compensates the multipolar effect. 
Taking a lattice-laser detuning of 5v = 0.11 MHz (white 
dotted line), light shift uncertainty |A//J'^(/, 
can be well less than 1 x 10“^® for an entire lattice 
laser intensity / shown in the plot. For 0.5% varia¬ 
tion in corresponding intensity range is reduced to 
/ < 12kW/cm^ by allowing fractional uncertainty to 
1 X 10“^®. The vector light shift arising from the el¬ 
liptical light may be canceled out by averaging the clock 
transitions with Zeeman substates m = ±1/2 24|. ^^^Yb 
and ^®®Hg with nuclear spin of 1/2 are the best candidate 
for this scheme because of the lack of tensor light shift. 

The sign of hyperpolarizability is not tunable for Sr. 
However, as shown by the region bound by a dashed 
rectangle in Fig. 4b, |Ar'f"'|/i/g'' < 1 x 10“^® holds for 
0 < / < 3.3kW/cm^ for a linear-polarized lattice (^ = 0) 
with 5v = 1.5 MHz, which makes Sr an attractive candi¬ 
date. Its very low Doppler temperature allows Sr lattice 
to operate at low intensity / ~ 2.3kW/cm^, where high 
merit factor Kgj. keeps the multipolar effect small and the 
hyperpolarizability effect does not come into play. Simi¬ 
lar low intensity optimization is applicable to Yb and Hg 
at the expense of lattice-trapped atoms or by applying 
deep laser-cooling on the — ®/b clock transitions. 

In case the hyperpolarizability is not tunable and/or 
experimental issues require to use linear-polarized lattice 
= 0), an optimal detuning 6vop may be determined 

= 0. Although the clock transi- 


as 


dAi'cjljSi'op) 


dl 


I=Io, 


tion may suffer from finite correction Ar'c(/opj the 

scheme would also work as the light shift becomes insen¬ 
sitive to lattice intensity A //1 for a certain range as seen 
in black-solid and dotted lines in Fig. |3l 


SUMMARY AND OUTLOOK 


We have proposed an operational magic frequency that 
makes light shift smaller than 1 x 10“^® for a wide range 
of lattice intensity as a result of cancellations of light 
shifts of different origin. Numerically demonstrated op¬ 
erational intensities with A/// > 40% offer a robust pro¬ 
tocol to guarantee superb reproducibility of optical lat- 
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tice clocks. It is noticeable that this intensity allowance 
can be larger than the intensity variation arising from 
the thermal motion of atoms in the transverse direc¬ 
tion A/// Ri A:BT’/-De(g), as it can be in the range of 
0.1 — 0.2 for atoms in thermal equilibrium 2^. Moreover 
all the tuning parameters, both the lattice-laser inten¬ 
sity / and the degree of circular polarization can be 
spectroscopically identified via the vibrational frequency 
n and the vector light shift, respectively, allowing one to 
define / and ^ in terms of frequency. In order to apply 
the proposed scheme, experimental investigations of the 
polarizabilities are crucial, where recently demonstrated 
clock reproducibility at 2 x 10“^® Q can be a power¬ 
ful tool. With precise determination of lattice intensity 
A///, e.g., by further cooling the transverse motion of 
atoms down to the photon-recoil temperatures [2l| , clock 
uncertainties of low 10“^® will be in scope. 
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